#include "parsing.h"
#include <math.h>
#include <float.h>

#define min(X, Y)  ((X) < (Y) ? (X) : (Y))
#define max(X, Y)  ((X) > (Y) ? (X) : (Y))


float hn(int n);
float sumc(int** r, int x, int* n);
float est_seq_chr(int r);
float fhat( int x,  int** a, float** c, int nind, float* var );
float theta_w(int** a, int* n , int ws,  int** r, int* s);
float theta_pi(int** r,int **a,float **c, int* n, int ws);
float sigma(int i, int j, int n);
float variance_tajima( int l, int nind, int* n, int** r, float theta );
